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ABSTRACT 

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts 
are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback 
times i > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200 < A < 2000 A) 
on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift 
approximately as Hubble parameter H~^{z). This is in accord with the theoretical expectation that 
the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial 
radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well- 
resolved sample of galaxies at z ~ 4 is 6/a = 0.65, suggesting that these Lyman break galaxies are 
not drawn from a spheroidal population. However the median concentration index of this sample is 
C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than 
to the values, C < 2 for local disk galaxies of type Sb and later. 

Subject headings: Galaxies: evolution — Galaxies: high-redshifts — Galaxies: structure — Cosmology: 
observations 



1. introduction 

An important goal of cosmology is to understand 
how galaxies evolve toward their current sizes and 
shapes. The basic framework of galaxy formation 
within the hierarchical Co ld-Dark-Matter (CDM) cos- 
mology was set out by iWhite fc ReesI l)1978t) . and 
has been refined by numerous N - body and semi- 
analytical studies fCole et al." '1994"; 'Kauffmann et al.' 
[W93; Somcrvillc & Primack 1999; White & Frcnk 1991). 
The formation of galactic disks within dark- matter ha- 
los was studied by .Fall fc Efstathioul 1)1980(1 . In this 
model, dark halos acquire their angular momenta via 
tidal torques, the angular momentum per unit mass of 
the baryons and the dark matter arc initially the same, 
and angular momentum is conserved as the baryons col- 
lapse and cool to form a disk. With these simplifying 
assumptions the baryons typically collapse by factors 
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of ~ 10 and the resulting disks have rotation curves, 
surface-density profiles, and scale radii similar to those 
observed. Further analytical studies have calculated 
the distribution function of dis k-galaxy sizes and the 
disk-galaxy size-redshift relation fP alcanton et al.(ll997t 
Mo ct al. 1998). Meanwhile cosmological N-body -I- hy- 
drodynamical simulations have formed disks that appear 
similar to spiral galaxies, but that tend to be too small to 
match present-da.y gala xies l|Navarro fc Steinmetdll99^ 
but see lEke et al.ll200'o|l . It is not yet known whether the 
source of this discrepancy lies in baryonic, dark-matter 
physics, or numerical issues in the simulations. 

If the general Fall & Efstathiou view of disk forma- 
tion is correct, there are several rather robust expecta- 
tions that are worth exploring through the observations 
of high-redshift galaxies. 

The size-redshift relation. The sizes of galactic disks 
forming at a redshift z should be a fixed fraction of the 
size of the dark-matter halo. The virial-radius of a dark- 
matter halo scales with redshift and virial velocity V^ir 
or virial mass Mvir as 



-,1/3 



100iJ2(^ 



lOH{z) 



(1) 



where 



H{z) = Hoprail + zf + + zf + Oa]1/2 (2) 

is th e Hubble parameter at redshift z ("Carroll et alJ 
Il992|) . Assuming that the exponential scale- length Rs 
of the baryonic disk scales with the virial radius, the 
sizes of disks are expected to scale with redshift as 
Rs oc H^^{z) at fixed circular velocity, or Rs cx H~^^^{z) 
at fixed mass. Observations are more suitable to track- 
ing evolution with fixed luminosity L, which may be ex- 
pected to fall somewhere in between these two functional 
forms. Previous observations have de monstrated that 
high-redshift galaxies are small (e.g. iLowenthal et alJ 
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Il997t iBouwens et al.ll2003ft . The larger sample of galax- 
ies provided by GOODS allows a closer look at the trend 
with redshift. 

The size distribution. If governed primarily by the an- 
gular momenta of their halos, the sizes of disk galax- 
ies should should show a log-normal distribution propor- 
tional to the dimensionless spin parameter: 
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where A is related to the total angular momentum J, en- 
ergy E, and mass M of the dark matter halo through 
the definition A = J\E\^^'^G-'^M~^/^. The predicted 
functional form for the distribution of disk-galaxy sizes 
appears to be a reasonable match to observations of 
present-day disk galaxies, albeit with a width cr that is 
signi ficantly smaller th an predicted by tidal-torque the- 
ory ijde Jong fc Lacevi r2000'l. It is sensible to imagine 
that a simple proportionality between galaxy radius and 
halo radius might hold even for elliptical galaxies, al- 
though models are less well developed in this case. It is 
less obvious that the size distribution should show the 
same functional form. 

The ellipticity distribution. The theory described 
above is specific to disk galaxies. It is thus interest- 
ing to test whether or not high-redshift galaxies are pre- 
dominantly disks. A relatively straightforward test is to 
compare the observed distribution of axial ratios to the 
distribution expected for disks or spheroids viewed at a 
random set of orientatio ns, ap plying the classic test out- 
hned bv lSandage et all l|197(T() . 

Throughout this paper we use a cosmology with pa- 
ramet ers fe, Ot.nt, ^^m i = 0.7,1.0,0.3,0.7, in the nota- 
tion of lCarroll et al.l l)1992t) . The sizes of galaxies quoted 
in this paper are in physical (proper), not co-rnoving , 
units. All magnitudes are on the AB system ijOkdlig?"^ . 

2. OBSERVATIONS AND MEASUREMENTS 

The Chandra Deep Field South (CDF-S) and Hub- 
ble Deep Field North (HDF-N) were observed with the 
Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) as part of the Great Observatories Ori- 
gins Deep Survey (GOODS). For the HDF-N, U-band 
observations were carried out using the Mosaic camera 
on the 4-m Mayall telescope. Details of the observa- 
tions, dat a reduction, and catalog generation are de- 
scribed in iGiavalisco et alJ l)2003air . Samples of galax- 
ies in broad bins of red shift were constructed using the 
[Giavalisco et al.' ('2003b') Lyman-break selection for z > 3 
and the Mobashcr et al. (2003) CDF-S photometric red- 
shifts (for two bins spanning the range z — 1 — 1.9 and 
z = 1.9 — 2.8). The z « 3 and 4 samples use data from 
both the HDF-N and the CDF-S. The z « 3 sample is 
drawn only from the HDF-N, where deep U-band images 
exist. Because the redshift selection is based entirely on 
photometry, there is undoubtedly some scatter of galax- 
ies from one redshift bin to another, and a few interlopers 
from lower redshifts in the high redshift bins. Based on 
simulations of template galaxies of various spectral types 
observed at the GOODS S/N, we expect the contamina- 
tion to affect only the tails of the observed size, elliptic- 
ity, and concentration-index distributions. Spectroscopic 
confirmations are of course highly desirable. 



Two different techniques were used to measure the 
radii of the galaxies. Half-light radii (denoted bv rh) wer e 
measured using SExtractor llBerti n fc ArnoutsI Il996j) . 
For this analysis SExtractor performs circular-aperture 
photometry using an analysis aperture radius that is 
ten times larger than the first radial-moment of the 
light distribution defined by connected pixels more than 
0.6 CT brighter t han the background in a smoothed de- 
tection image. iPetrosia n ( 1976!) radii were measured 
using the prescription adopted for the SDSS project 
ijStoughton et al.ll2002|) . The ratio of the semi-minor to 
semi-major axis b/ais measured by SExtractor using sec- 
ond moments of the light distribution within the analysis 
aperture. 

To compare the sizes of galaxies in different redshift 
intervals, we confine our samples to non-evolving rest- 
frame luminosities between 0.7L* and 5L*, where L* is 
the characte ristic lumin osity of a z = 3 Lyman break 
galaxy from iSteidel et a l. (1999). Five samples in dif- 
ferent bins of redshift (z w'l.4, 2.3, 3.0, 4.0, and 5.0) are 
constructed. For each sample we select galaxies within 
the desired luminosity range based on their apparent z- 
band magnitudes, after accounting for cosmological ef- 
fects and fc-corrections assuming a typical Lyman-break 
galaxy spectral-energy distribution (SED). The adopted 
SED is from the 1999 version of the Bruzual & Charlo^ 
(^993) solar-metallicity models, with a Salpeter initial- 
mass function, constant star-formation rate, and age of 
144 Myr This is attenuated assuming E{B ~ V) = 0.14 
with the lCalzetti et all 1)20001) extinction curve to match 
the mean UV spectral s lope seen in the z ~ 3 sam- 
ples (jAdelberger fc Steid cl 2000). The luminosity range 
translates to a z-band magnitude range of 23.1 to 25.3 
for the z ~ 4 sample, with that magnitude range shift- 
ing by -2.41, -1.12, -0.64, and -1-0.46 magnitudes for the 
samples in the redshift bins at z ~ 1.4,2.3,3.0, and 5.0, 
respectively. After rejecting sources classified as stellar 
(SExtractor CLASS_STAR> 0.9) or with unrehable pho- 
tometry (SExtractor FLAG> 3), the sample sizes (for the 
lowest to highest redshift bins, respectively) are 21, 73, 
140, 386, and 153 galaxies. 

3. RESULTS 

The galaxy sizes are measured in the band closest to 
rest-frame 1500A. At z = 3,4, and 5, the central wave- 
lengths of the F606W, F775W, and F850LP band fihers 
fall at rest-frame A = 1500 ±40 A. For the lower-redshift 
samples the size measurement was made in the F435W 
images corresponding to A ~ I8OOA at z ~ 1.4 and 
A ~ I3OOA at z = 2.3. The sizes for any of these samples 
do not change appreciably when measured through the 
next redder bandpasses. 

The measured size-redshift relation for our samples is 
shown in Fig. ^ The observed points are uncorrected for 
biases or incompleteness (see 2J j ^-nd the errorbars rep- 
resent the standard error of the mean for each sample, 
ignoring possible effects of clustering (jSomerville et alJ 
i2003.') or crowding. The measurements are roughly con- 
sistent with constant angular sizes above z = 2.5. For 
comparison, in Fig. ^ we show three size-redshift rela- 
tions. The blue line is the angular-size redshift relation 
for a standard measuring-rod in the WMAP cosmology 
(jSpergel et al.ir2003 ). In this cosmology, angular diam- 
eters increase by 23% from z = 3 to z = 5 (for objects 
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with fixed proper diameters) . The two hnes that decrease 
with redshift are the expected scahngs if size at fixed lu- 
minosity tracks size at fixed virial velocity r cx H~^{z) 
or fixed mass r cx H~^/^{z). The mean half-light radius 
of 0.24" at z « 4 corresponds to 1.7 kpc. 



for flattened systems is indeed robust. Visual inspection 
of the images reveals very few galaxies that have mor- 
phologies reminiscent of nearby disk galaxies. The ob- 
served flattening distribution is probably a reflection of 
the tendency for many of the Lyman-break galaxies to 
have a few concentrations of light of nearly equal bright- 
ness rather than a central dominant concentration. 
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Fig. 1. — Size vs. redshift relation. Mean SExtractor half-light 
radii are plotted with errorbars indicating the standard error of the 
mean (i.e. the sample standard devation divided by the square-root 
of the sample size). The solid blue curve shows the expected trend 
in the WMAP cosmology if physical (proper) sizes do not evolve. 
The dashed red curve shows the trend if sizes evolve as H~^{z) 
and the dotted green curve shows /f-2/3(z). The curves are all 
normalized to the mean size at z si 4 (approximately = 1.7 
kpc). 

The size distributions of galaxies in the different red- 
shift bins are shown in Fig. |21 The curves show a "no- 
evolution" model tuned to match the observed size dis- 
tribution at z = 4. This model accounts for selection 
and measurement biases (including cosmological surface- 
brightness dimming) as described in 21 The galaxies 
in the two low-redshift bins are clearly systematically 
larger than the no-evolution model would predict. How- 
ever, given the sample sizes and measurement biases, the 
trend between redshifts z = 3 and z = 5 is not easily dis- 
tinguished from the no-evolution model. 

Figure 12^ compares the observed distribution of cUip- 
ticities for the B-dropout z « 4 sample to simulations. 
The observed population has a larger proportion of galax- 
ies with axial ratios 6/a < 0.4 than expected for a purely 
spheroidal population. A Kolmogorov-Smirnov test indi- 
cates a likelihood of 57% that observed distribution could 
have been drawn from the model disk distribution, while 
the probability is 5 x 10~^ that the observed distribution 
could have been drawn from the model spheroidal dis- 
tribution. Similar results (at lower significance) are ob- 
tained for the z ^ 3 and z ~ 5 samples. However, these 
statistical probabilities are sensitive to the size and mag- 
nitude distribution of galaxies fed into the simulations. 
While we have attempted to match the observed distri- 
bution of magnitudes and sizes (and colors) 'At z — A (see 
Fig. 2), more work is needed to verify that the preference 
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Fig. 2. — Observed size distributions. Data for the differ- 
ent samples are shown with Poisson error bars. The curves show 
a "no-evolution" model of the size distribution, tuned to match 
approximately the observed distribution at z fa 4. The model 
has a log-normal distribution of radii with a peak at 2.1 kpc (for 
L = 2.9L*) and a width cr = 0.3. We adopt the size— luminosity 
relati on (i? oc L^/^), observed for local disks (dc Jong & Lace3 
120001) . The models include a 50/50 mix of oblate ellipticals and 
flat exponential disks. Model galaxies have been inserted into the 
images in a Monte-Carlo fashion and detected and measured using 
SExtractor in the same way as the real galaxies. Therefore the 
model curves faithfully represent the effects of both sample incom- 
pleteness and measurement biases. The input distribution (prior 
to selection and measurement biases) for the z fa 4 bin is shown as 
the dashed curve. Objects classified by SExtractor neural-network 
as stellar at greater than 90% probability are excluded from the 
samples (both the data and the simulations). In practice only a 
few objects with < 0.1" are removed. 



The concentration-index distribution for the same sam- 
ple is shown in Fig. ^p. The concentration index here 
is defined as 5 log(r8o/?'2o) where rgo and r2o are radii 
containing 80% and 20% of the galaxy flux, respectively 
ijConse licc 2003). The observed galaxies span a broader 
range than the ideal disks and spheroids that populate 
the simulations. The z « 4 population tends to be more 
centrally concentrated than pure exponentials. Nearby 
spiral galaxies tend to have central concentrations of light 
as well, but these also tend to be red, and virtually invis- 
ible in the rest-frame ultraviolet where the current mea- 
surements are made. 
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Fig. 3. — (a) EUipticity distribution for 2 4 galaxies com- 
pared to simulations for pure disk and pure spheroid populations 
(the blue-dashed and red-solid curves, respectively). For the simu- 
lations, the disk sample (viewed from arbitrary directions) is drawn 
from a population of oblate optically-thin spheroids with a Gaus- 
sian distribution of intrinsic axial ratios with mean b/a = 0.05 and 
(T = 0.01. The spheroids are drawn from a population of oblate 
spheroids with intrinsic axial ratios uniformly populating the range 
0.3 < b/a < 0.9. The samples for both observations and simula- 
tions are limited to galaxies with 23 < ^775 < 25 and SExtractor 
half-light radii > 0.2". 

(b) Concentration-index distribution for for z fa 4 galaxies com- 
pared to simulations for pure disk and pure spheroid populations 
(the blue-dashed and red-solid curves, respectively). Samples are 
the same as in (a). 



4. SELECTION EFFECTS AND BIASES 

In this section we highlight the issues of bias and in- 
completeness and the tests that have been carried out to 
verify the results of the previous sections. 

Incompleteness. The completeness limits of the 
GOODS images have been estimated by inserting artifi- 
cial galaxies with a wide range of half-light radii and mag- 
nitudes into the GOODS z-band image, and re-running 
SExtractor. The simulated galaxies have either r^/^-\iew 
or exponential surface-brightness profiles, with the ax- 
ial ratio distributions described in the caption to Fig. 
3. The completeness li mits as a function of size and 
magnitude are shown bv iGiavalisco et all l)2003a(l . At 
-^850 = 24.7, corresponding to the mean magnitude of 
the z « 4 sample considered in this letter, the sam- 
ple appears to be 80% complete for galaxies with half- 
light radii rh = 0.5" and 50% complete for galaxies with 
Th = 0.9". The z > 3 samples are affected by incom- 
pleteness for galaxies with sizes ^ 3.5 kpc; for the 
lower-redshift samples there is no equivalent bias against 
detecting galaxies with small sizes. Therefore the most 
robust result is the observation that analogs of the galax- 
ies with small half-light radii in the z > 3 bins are much 
less common at the lower redshifts. Without careful 



modeling (or deeper images) it is not possible to rule 
out the possibility that galaxies with larger sizes exist 
at the higher redshifts. It is worth noting that a pure 
exponential disk galaxy with a luminosity L* and a scale 
length of > 5 kpc would be unlikely to be detected at 
z > 3 in the GOODS z-band images. 

Radius-measurement biases. The results shown in Fig. 
1 and 2 are for SExtractor half-light radii. The trends 
and distributions look virtually identical when Petrosian 
radii are used. Based on the simulations it is clear that 
both radius measurements become biased at faint mag- 
nitudes and small sizes. For example, for galaxies in the 
magnitude range of our z « 4 sample, a galaxy with an 
intrinsic r/j — 0.5" would have a measured = 0.36". 
Measurement biases will thus pull in the tail of the size 
distribution; the true size distributions at z > 3 are al- 
most certainly broader and peaked at slightly larger sizes 
than the observed distributions shown in Fig. 2. The 
model distribution for the z ~ 4 sample provides an es- 
timate of the importance of this bias. The true distribu- 
tion of half-light radii input to the simulation is shown 
as the dashed line. For this input distribution the effect 
of selection and the size-measurement bias are minor. 
A distribution with a more significant tail to large sizes 
would show more of an effect. We have tried broader dis- 
tributions and distributions with a larger mean size and 
found them to predict more large galaxies than observed 
in the z « 4 sample. 

5. SUMMARY AND CONCLUSIONS 

Samples of galaxies within a fixed and fairly narrow 
range of rest-frame UV luminosities have been compared 
from redshifts z 1 to 2: 5. The sizes of high-redshift 
galaxies appear to evolve roughly as H~^{z) in agreement 
with the rather robust expectation from hierarchical the- 
ory. This general conclusion is subject to the caveat that 
we may be seeing different kinds of galaxies at the dif- 
ferent redshifts. Our sample selection necessarily targets 
UV-bright objects; further studies on samples selected in 
the near-infrared or in subsets of morphology are clearly 
warranted. For the galaxy sample discussed here, the size 
distribution at z ~ 4 is reasonably well represented by 
a log-normal distribution, but with a smaller mean size 
than observed for samples of nearby luminous galaxies. 
Measurements of the ellipticity and concentration-index 
distributions suggest that Lyman-break galaxies at z ~ 4 
represent a mix of morphologies, with some tendency to- 
ward flattened systems. 
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